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, We theoretically and numerically studied ordering of "colloidal ionic clusters" on periodic sub- 

■ strate potentials as those generated by optical trapping. Each cluster consists of three charged 

Lj^ ' spherical colloids: two negatively and one positively charged. The substrate is a square or rectan- 

gular array of traps, each confining one such cluster. By varying the lattice constant from large to 
small, the observed clusters are first rod-like and form ferro- and antiferro-like phases, then they 
bend into a banana-like shape and finally condense into a percolated structure. Remarkably, in a 
broad parameter range between single-cluster and percolated structures, we have found stable super- 
complexes composed of six colloids forming grape-like or rocket-like structures. We investigated the 
C/j ■ possibility of macroscopic pattern switching by applying external electrical fields. 



oo 



- 1—1 

X 

5-H 



I. INTRODUCTION 
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It is generally considered that the study of colloidal systems started with the work of Thomas Graham, around 
1860, although substances such as the "purple of Cassius" , which turns out to be colloidal gold, were used by glass- 
makers and porcelain manufacturers as early as the XVIIth century [l|. For a long period, those systems were 
considered as "unmanageable" (see e.g. the preface of Ref. [U), due in particular to the difficulty to characterize 
■ them precisely (size, shape, charge etc) and also to the lack of theoretical tools to understand effective interactions. 
| From this perspective, a decisive step forward occurred in the 1940s with the work of Dcrjaguin, Landau, Verwey and 
Overbeek 2| > which helped turn the initial empiricism into a more solid body of knowledge. More recently, progress 
in experimental manipulation techniques promoted colloids as interesting model systems to investigate a large gamut 
of fundamental physical phenomena: phase ordering, nucleation, effects of confinement and reduced dimensionality, 
interfacial phenomena, and various non-equilibrium problems including glass formation, to name but a few. Indeed, 
compared to atomic systems, colloids offer four advantages : first, they are characterized by small elastic constants 
and therefore they are easy to manipulate ( "soft matter" ) ; second, they are more easily visualized since their typical 
length scale is within the range of visible light; third, the corresponding time scale is no longer in the pico-second 
range, but of order Is; fourth, their effective interaction potential may be easily tuned to some extent by changing 
^ \ some control parameter such as the concentration of a given solute Q . Such features stem, directly or indirectly, from 
the large size (on an atomic scale) of the mesoscopic constituents. Hence, colloids can be considered as "big atoms" 
but with controllable interactions, and it is noteworthy that they may furthermore exhibit new structures, that 
d . do not seem to have any atomic analogue Q. 

Of particular interest here is a class of charged colloidal spheres in the presence of optical trap arrays, where a wide 
variety of crystalline states has been reported f6Hl"2l]. On a two-dimensional modulated substrate, multiple charged 
colloids can be confined in sufficiently strong traps, and thereby form ?i-mers which exhibit an orientational degree of 
freedom. It has been established from experimental, numerical and theoretical approaches [6l-ll2| that such systems 
may show remarkably rich orientational ordering, and lead to the formation of so called colloidal molecular crystals. 

Previous approaches considered situations where the colloids forming trapped n-mers are like-charged objects. On 
the other hand, the 'molecules' under study here are made up of oppositely charged colloids, with equal charges 
in absolute value. The focus will be on a square lattice of traps, including also distorted geometries where all 
distances along one of the principal axis of the square is scaled by a given factor a, thereby producing a rectangular 
lattice. We will restrict our attention to the ground state, leaving thermal effects for further studies [Hj]. Under 
these circumstances, it has been shown that in the dipolar case (with trapped clusters each formed of two oppositely 
charged colloids) the resulting orientational behaviour is somehow trivial: the observed phase is made up of ABAB- 
type stripes, where the dipoles in stripe A are aligned along the stripe (up), while those in stripe B have opposite 
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orientation [12j. Furthermore, the stripe structure cannot be tuned by modifying the lattice constant or salt content 
in the solution. Here we consequently consider clusters with three colloids, two negatively and one positively charged, 
which are more amenable to an external control, and display a wealth of ordering patterns, which is the center of our 
interest. The paper is organized as follows. The model and the methods are defined in section|TTl where the differences 
in modelling with previous works are outlined. Behaviour on square and rectangular lattices with isotropic parabolic 
confinement potential is studied in section IIIII and in section IIV1 we discuss the relevance of these results in the case 
of a more realistic cosine trap potential. All these results pertain to equilibrium configurations, while orientational 
pattern modifications induced by an external electric field is finally addressed in section [V] In the conclusion, we 
summarize our main findings and discuss perspectives and open problems left. 



II. MODEL AND METHODS 



We consider a mixture of negatively and positively spherical and homogeneously charged colloids, with stoichiometry 
2:1 (twice more negative than positive macroions, which is equivalent to the situation where all charges are reversed). 
Those mesoscopic objects bear a total charge that is neutralized by an ensemble of micro-ions. We assume that 
counter-ions and electrolyte micro-ions do not differ. The colloids are confined in a two dimensional plane by the 
action of laser beams. This plane is furthermore corrugated by additional preferential confinement: the colloids, 
subject to gradient forces and light pressure [H, [l4[ tend to gather in the regions of highest laser intensity. It is then 
experimentally possible to create a 2D periodic substrate of traps, with variable geometry. It should be emphasized 
that the micro-ions are not sensitive to confinement, so that our two dimensional colloidal system is immersed in a 
standard electrolyte, that mediates between the charged colloids a screened Coulomb interaction of the form [1, [H, EH] 

^ = A'£ S «i=f. (1) 

The energy scale K, which contains a geometric prefactor together with the squared effective charge of the colloids 
[17) . is not relevant for our purposes, since we are interested in ground state configurations. In Eq. (fTJ), the sum 
runs over all colloidal pairs with interaction sign s,j = ±1 distinguishing like- from oppositely charged colloids, and 
inter-particle distance between macroions i and j; with I we denote the intertrap distance, see Fig. [TJ We therefore 
implicitly assumed that oppositely charged colloids at contact also interact through a screened Coulomb form, whereas 
such a functional dependence is better suited to describe the far-field behaviour [18| . We note however that as long as 
we are comparing phases with an equal number of doublet +/- at contact, such configurations contribute a constant 
to the total energy, so that considering more exact forms of interaction potential at short scale would not alter our 
conclusions. In Eq. ([I]), 1/n denotes the Debye length, that can be tuned changing the electrolyte density. Strictly 
speaking, Eq. (JXJ) is not an energy but a free energy, since it includes an entropic contribution from the micro-ions 

.19]- 

In addition to inter-particle forces (two-body, and pairwise additive), the colloids feel an external (one-body) 
potential due to laser modulation, and it is the subtle interplay between both kinds of forces that selects the type of 
ionic complexes that is formed in each trap, and their relative arrangement from trap to trap. In [HQ, it was assumed 
that the colloidal complexes in a given trap (dimers, trimers etc) form a rigid object, with a single orientational degree 
of freedom, which effectively allows to omit the laser-induced potential in the analysis. The situation with oppositely 
charged colloids however differs: take an isolated triplet -/+/-, in the aligned geometry which minimizes its self energy. 
Such an object invariably bends under strong confinement thereby forming a 'banana-like' object. The shape of the 
banana depends on the parameters (salinity, inter-trap distance etc) , and should be determined self-consistently. We 
therefore do not lump substrate potential effects into a unique quantity as in 0, Q, but explicitly consider each 
colloidal degree of freedom in the analysis. This rules out the discrete angle approach invoked in [9J, where only a 
restricted set of pos sible orientations were considered. We also note that already at the level of like-charged colloids, 
our recent work [10( evidenced some subtle shortcomings of the rigid n-mer approach. 

The system obviously bears some resemblance to the molecules on the atomic scale. However, in our model the 
binding is the result of the attractive Yukawa interactions and the confinement potentials and is therefore relatively 
weak compared to the strong covalent bonds typical in real molecules. We have consequently chosen to call the n-meric 
colloidal structures "clusters" or "complexes" rather than "molecules" as has been done elsewhere (see 0413 )• ^ n 
spirit, our study then belongs to a more general line of research pertaining to colloidal clusters [2(j. The previous 
references though generally deal with large assemblies, whereas our clusters are made up of a small number of colloids. 

In all our considerations, the ratio of the total number of colloids to the total number of traps is exactly three, 
therefore under strong enough confinement, there are exactly 3 colloids per trap and the ground state is paved with 
trimers ( — I — ) of varying shapes and orientations. These mesoscopic objects bear a net charge, which is neutralized 
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by micro-ions in solution. It has been realized [g, [10( that it is not necessary to invoke local anisotropies of the 
confinement potential to account for the observed ground states: there is an obvious source of anisotropy in the 
lattice geometry itself; however, in the vicinity of a trap minimum, the confined particles feel an isotropic (2D) 
potential. In this case, the screened Coulomb potential alone is able to select preferred orientations for the clusters 
(whereas as alluded to above, the effective spin model considered in [9( restricted possible orientations to the principal 
directions of the light lattice). Here we wish to avoid introducing an orientational bias through the trap potential, 
and therefore initially consider the parabolic confinement. The total dimensionless energy of colloid i reads 



where Si denotes the distance between colloid i center of mass and the center of the trap to which it belongs, and 
A measures the relative importance of parabolic confinement over Coulombic forces. Each triplet of colloids belongs 
to a 'native' trap, so the colloids are not allowed to hop from one trap to another. Such a limitation precludes the 
sliding states studied in [2l[ under the influence of an external electric field. We note that the experimental situation 
is somehow intermediate between the point of view ado pte d here - a priori mostly relevant for strong confinement 
conditions -, and that used in Q, see the discussion in jlOl]. For completeness, we have also performed simulations 
with a 'realistic' potential in the spirit of [7| to assess the validity of the isotropic approach. 



FIG. 1: (Color online) Sketch of two clusters - each made up of three colloids - in two different potential traps, with definition 
of the relevant quantities. The center of a blue area denotes the confining potential minimum and I is therefore the intertrap 
distance. In all the paper, the same color code will be used: black (resp. white) dots for positively (resp. negatively) charged 
colloids. 

At constant temperature, three parameters govern the behaviour of the system: the relative magnitude of the 
confinement with respect to the Yukawa interaction A, the ratio of the colloid radius to the Debye screening length 
koI, and the ratio of the lattice spacing to Debye length kI. In the rectangular case, the lattice aspect ratio a should 
be included. A trimer in two dimensions has six degrees of freedom, but only four of them are independent assuming 
the colloids are in contact. The arrangement of colloids in the vicinity of a trap center can therefore be described 
by four quantities (see Fig. [IJ: the position of the cluster's center of mass relative to the trap minimum is described 
by r and 7, the relative orientation of its axis with respect to a fixed direction (say a line joining the centers of two 
neighboring traps) by the angle 1? and the amount of bending by (j) {<j) = it for rod- like and 4> < tt for bent structures). 

Ground state configurations were obtained by energy minimization (EM) technique and tested by standard Monte 
Carlo (MC) annealing simulations. In order to reduce the number of minimizing parameters in EM, an assumption 
on the ground state structure was preset. For a given rectangular lattice we decided to partition the lattice to 2 (A-B) 
or 4 (A-B-C-D) sub-lattices, which in general form a checkerboard or a stripe (chain) pattern. Since each colloidal 
cluster is parametrized by four quantities (see Fig. [TJ, a p-partite (p = 2,4) assumption in EM would correspond 
to 4p minimizing parameters. We used standard numerical routines (simplex and quasi-Newton algorithm (22j) for 
finding the (local) minimum; by repeating the EM procedure several (~ 10 ) times over 4p different initial estimates, 
we eventually converged to a global minimum. However, in extreme parameter regimes (i.e., 1/ d^> 1 or close to phase 
transitions), the convergence became slower due to many almost degenerate local minima and the ground state was 
identified among selected candidates - by directly comparing their energies. 

The 2- and 4-partite constraint used in EM was further assessed by standard MC annealing simulations. We 
checked all interesting regions of the phase diagram by running several MC simulations starting from different initial 
configurations (both minimal or close-to-minimal EM structures and random configurations) and used the annealing 
method with exponentially decaying temperature profile. We tried several system sizes with 2 x 2, 4 x 4 and 16 x 16 
number of traps using periodic boundary conditions. As a result, we identified and confirmed all shapes obtained 
in EM. However, due to the tiny energy differences and large energy barriers between local minima, long relaxation 
times (i.e., 10 8 MC steps for 4 x 4 system array) were needed to obtain exactly the same orientational orders as in 
EM. Hence, the system easily gets trapped at low annealing temperature in some local metastable phase composed 
of domains of various ground state patterns (23|. 
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III. SQUARE AND RECTANGULAR LATTICE CONFINEMENT 

Before embarking on a detailed study, we first give an overview (Fig. [2]) of the various composite objects that we 
have observed in our simulations. When the traps were well separated, we observed isolated trimers, either straight or 
banana-shaped B. For very small trap separation I, percolated chain-like structures C rather than clusters are stable. 
In the intermediate region "super-clusters" composed of six colloids appear. They are formed by two neighbouring 
trimers forming the super-cluster by minimizing their total energy (paying the penalty due to the confinement while 
reducing the electrostatic interaction energy) [24|. We observed two typical shapes and coined them grape Q and 
rocket 1Z. In Fig. [3] we further provide a summary of different orientational orderings and propose a nomenclature to 
classify them. 
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FIG. 2: Summary of the most frequent composite clusters encountered in our study (from left to right): straight trimer, banana 
B, grape Q, rocket 1Z. 
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FIG. 3: Sketch of the different ordering patterns that will be reported (from left to right): anti-ferro A, chain C, ferro T, tilted 
ferro T* and ladder L. Dashed lines represent a square (or rectangular) lattice of traps. 



A. Square lattices 



Figures S] and [S] show the type of orderings which occur on a square lattice (a = 1), at two different pinning 
amplitudes (embodied in the parameter A). On panel a) we show the phase diagram separating the regions of 
stability of phases and the corresponding snapshots of the typical cluster arrangements. The panels b) and c) display 
the variation of three quantities characterizing the cluster shape and their orientational ordering. The transitions from 
straight to bent shapes are clearly seen on bifurcating red curves. It should be stressed that </> and 2ir — cf> represent 
the same shape and the two branches do not indicate two different structures. The black curves showing the cluster 
orientation have two branches with -di and $2- The special case is the A (antiferro) phase, where ($i,$2) = (0, 7r/2). 
For straight clusters the angles $ and 2ir — are equivalent, but this is no longer true if they are bent. Therefore, 
the region with four distinct branches on Fig. [5] represents a 4-partite structure with i?i...i?4. The values of r (green 
curves) tell us whether the clusters are displaced from the trap minima (r ^ 0) or not (r = 0). The values of the 
fourth defining quantity 7 proved rather uninformative and were not included in the figures for clarity reasons. 

At the stronger confinement (Fig. 2]) we found three phases: isolated triplets in the anti-ferro like arrangement A, 
bent (banana) triplets B whose angular arrangement varies continuously to anti-ferro A upon increasing l/d, and a 
percolated chain pattern (C) at very low l/d and small screening (lower left corner of Fig. 0Ji)). At weaker confinement 
(Fig. [5]), two additional phases appear: rocket 1Z and grape Q, which are both non-percolating, and made up of the 
repetition of 2-triplet (6 colloids) objects. 

In Fig. 21 all ground state configurations have a 2-partite lattice structure, which has also been confirmed by the 
unconstrained MC annealing simulations. In the limit of small macro-ions, l/d 3> 1, we find stable isolated triplet 
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FIG. 4: (Color online) a) Ground state phase diagram at a relatively strong substrate potential (A = 5) for colloidal trimers 
on a square lattice. Typical snapshot configurations are calculated using (l/d,K,l) = (7,2) for the B phase, (15,4) for the A 
structure, and (5, 2) for the C phase. The lines between the different phases are shown as a guide to the eye. b) Values of the 
three quantities characterizing the cluster shape: (red dashed (upper): r, black solid: •d and green dot-dashed (lower): ip), as 
a function of the ratio l/d at constant nl — 7. c) Values of the same three quantities as a function of kI at constant l/d = 8. 



configurations with quasi degenerate orientational preferences. For example, the anti-ferro A phase is slightly more 
favorable compared to the tilted ferro J 7 * phase: at the (l/d, nl) = (15, 4) point in the phase diagram the corresponding 
energies are E A /K= -2.4250 and E?* /K = -2.4247. 

In Fig. [5] there are regions in the phase diagram where the 4-partite configurations appear to be favored over 
bipartite. The energy difference between 2- and 4-partite ordering, however, is extremely small. For example, in 
the representative rocket 7Z and grape Q phases shown in Fig. [Sfa) the energies are respectively E 4p /K = —0.5888, 
E 2p /K = -0.5870 and E 4p /K = -2.723145, E 2p /K = -2.723143. Such small differences, E 4p - E 2p , imply that 
these ground state configurations become unstable at small but finite temperatures and the formation of domains is 
expected in experiments and simulations. Similar effect is expected also for the A ground state, which should interfere 
with nearest (in energy) T* phase, as in the previous A = 5 case in Fig. [5] On the other hand, the grape Q phase in 
the l/d 3> 1 limit is well separated from the competing A phase (Ej^/K — —2.645310 for the same nl and nd values). 



1. Antiferro-like phase 

The occurrence of the A phase on the right hand-side of Figs. 0] and [5] can be rationalized following similar lines 
as in Ref. Q. For large enough values of kI and even larger l/d, it is legitimate to restrict the summation in Eq. ([2]) 
to nearest neighbor traps (see (Tol j for a more complete discussion). It is also legitimate to assume that triplets align 
in straight configurations. The question is therefore to understand the relative orientation of triplets, from trap to 
trap. Elaborating on the remark, put forward in [25|,[26j], that a non isotropic charge distribution in an electrolyte, 
creates an electric potential that is anisotropic at all scales and that does not have any simple 'multipolar' symmetry, 
we write the angular dependence of the far field potential of a triplet as [8j 

Vtr = 2cosh(2KiisinV') — 1 (3) 

where the direction ip is defined in Fig. [SJa). Consequently, the relevant angular dependent triplet-triplet interaction 
potential may be written (see Fig. [6]-b)): 



V tr /tr = [2cosh(2K<isin?/!i) — 1] [2 cosh(2Kdsin-02) — 1] 



(4) 
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FIG. 5: (Color online) Same as Fig. [4]but for a relatively weak substrate potential (A = 2). a) Typical snapshot configurations 
are calculated using (l/d,nl) = (7,8) for B, (7,4) for TZ, (5,2) for C, (20,8) for G, and (15,2) for A phase. 




FIG. 6: Definition of the notations used to compute triplet-triplet interaction in the (aligned) straight case. 



On a square lattice, and assuming bipartite structure (which is backed up by simulations), we then have to consider 
triplet-triplet interactions along both principal axis of the lattice. This leads to the following energy 

[2 cosh(2K<isin?/>i) — 1] [2 cosh(2Kdsmip2) — 1] 
+ [2 cosh(2ft<icos?/>i) — 1] [2 cosh(2«;<icos ^2) — 1] • (5) 

The minimum is reached for (^1,^2) = (0,tt/2), irrespective of nd, which exactly coincides with the A phase con- 
figuration. The situation changes on a rectangular lattice, see below. On the other hand, the A — > B transition 
line observed in Fig. [4] can be explained as a one-cluster effect: in the absence of confinement, the trimers take 
on the straight shape due to the repulsion between the two positively charged colloids. However, by increasing the 
confinement or because of enhanced screening, the relative importance of the electrostatic repulsion decreases. We 
expect that for A and d fixed, there exists a critical inverse screening length n* beyond which the straight trimer 
bends (and ultimately forms an equilateral triangle if k is very large). Since it is a one-triplet instability, k* should 
be independent of lattice spacing I which translates into a straight line with slope K*d in a (nl, l/d) plot such as Fig. 
2Ja). This corroborates our numerical finding with a straight separatrix between A and B phases. For A = 2, the 
argument does not hold since no transition from straight to bent shapes takes place. 



7 



2. Grape phase 

The stability of the grape phase Q can be rationalized by comparing the self-energy per colloid 

_ d A e- Kr ^ 

E S elf - g 2j S V r . . 



(6) 



of different composite complexes in the absence of any trap potential. In Fig. [7] we have compared the self-energies of 
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FIG. 7: (Color online) Coulombic self-energy of 6 colloids (2 positively, and 4 negatively charged) forming various complexes. 

the following six-colloid objects: grape, rocket (with various angles e, as defined on Fig. [2J and two distant trimers. 
Among these, the two distant trimers are the most stable in the bare Coulomb limit (n — > 0) with E se if — > —1/4. 
However, as soon as Kd > 2.2 • 10~ 2 , a fairly small number, the grape geometry is found to be the optimal one. 
Therefore, in almost the entire range of relevant values, the grape is energetically the most stable isolated object 
consisting of six colloids. Of course, on a corrugated substrate, the confining potential influences the above scenario 
by penalizing extended structures, and the inter-cluster interaction affect the picture as well, however, there is still a 
substantial region on Fig. [5] and [8] (see below) where the grape phase remains stable. 



B. Rectangular lattice of traps 

Distorting lattice geometry by changing the aspect ratio away from the square shape (a = l y /l x < 1) induces new 
orderings to the list of ground state configurations reported above. Figure [8] shows ground state phase diagrams for 
a = 0.8 and two pinning potentials, where three new phases emerge with respect to the square (a — 1) case: a ferro 
T phase and a tilted ferro J-* phase, which evolves into a familiar banana B phase at smaller l/d and A = 5, and 
a percolating ladder C phase in a weaker confinement (A = 2). The latter one consists of complex subunits where 
each trimer touches with three (out of four) neighboring trimers. Both pinning regimes involve also a distorted (or 
bent) rocket TZ phases, which we identify as a variant of the existing (a — 1) phase rather than a new phase itself, 
and a grape Q phase with 2- and 4-partite structures. Actually, in Fig. EJb), almost the whole parameter space is 
minimized by a 4-partite grape Q configuration (27l |. 

On a rectangular lattice, we can follow a similar line of reasoning as above in section fill A| where the occurrence of 
the A phase was justified. In contrast to the situation on the square lattice, on the rectangular lattice the number of 
nearest neighbor traps decreases to two and in order to assess the stability of the T phase, it is sufficient to consider 
a single triplet-triplet interaction of the form Eq. (|4]). In this case, the lowest energy configuration corresponds to 
parallel triplets (tpi — ^p2 — 0) forming a ferromagnetic-like arrangement along the longest direction of the lattice 
unit cell. This is indeed what is observed in Fig. [5Ja) (see the T region on the right hand side), but not in Fig. 
[8jb), because in this latter case, confinement strength is too weak to prevent neighboring traps from sharing their 
macroions. We further note that the T* — > B transition visible in Fig. |S][a) falls in the same category as the A^r B 
seen in Fig. 2) it is a one-cluster instability due to the competition between Coulomb and confinement forces. It is 
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FIG. 8: Ground state phase diagram for colloidal trimers on a rectangular lattice with a = 0.8 and a) A = 5 or b) A — 2. 
Typical snapshot configurations in a) are calculated using (l/d, kI) = (10, 8) for B, (7, 4) for 1Z, (7, 2) for Q, (15, 6) for T*, and 
(20, 2) for T phase; in b) (l/d, nl) = (7, 8) for 11, (5, 6) for £, (7, 2) for 2-partite Q, and (10, 2) for 4-partite Q phase. 



hence not surprising to find that the slopes of both separatrices are close: they are quantified by the same critical 
screening coefficient n* introduced above. 



IV. TOWARDS A MORE REALISTIC TRAP POTENTIAL 



On a realistic trap, when allowing the particles to hop, we can expect to significantly extend the stability of the 
grapes, at least as soon as nd exceeds a small threshold. This is supported by the evaluation of the self energies on 
Fig. [7] and we further analyzed the stability of the grape phase by explicitly introducing a cosine potential, which 
provides a reasonable approximation for a realistic optical confinement [2j: 

yreal = _ A_ ( cos ( 27ra //) + C0S (2TTy / Ctl)) . (7) 

The prefactor 2tt 2 is chosen such that both realistic and simplified potentials behave in the same way, as (S/l) 2 — 
(x 2 + y 2 )/l 2 close to the trap minimum, when one chooses A = A r . 

In Figs. |U[b,c,d) we show ground state configurations for various pinning strengths A r using the same parameter 
values (l/d,nl,a) as in Fig. [9ja) where we show the four most representative parabolic-trap configurations. We 
choose A r such that both potentials become similar in amplitude at the point of maximum Vf eal (see left panel of 
Fig. [9]). As A r gets smaller (going from Fig. [9jb) to Fig. [9jd)), the optimal configurations indeed evolve to the 
grape Q geometry. In the limit A r — > A, i.e., when the potentials share the same behaviour close to the trap center, 
the lowest configurations are always grape Q phases. Although having the same internal structure, these grapes are 
formed in realistic trap minima rather than at trap midpoints. In addition, the orientation of grapes is selected by 
the symmetry of the underlying realistic potential Vf eal , which is lowered with respect to the isotropic parabolic trap 
(the inset in the left panel of Fig. [HI shows the equipotential lines of Vf eal ). In fact, the realistic trap square symmetry 
is manifested in all configurations shown in Figs. [9][b,c,d). 

From Fig. \§\ we see that switching to realistic trap potential (Vf eal ) preserves some typical cluster shapes (e.g., 
straight trimer, grape Q, rocket 1Z) from the parabolic case, while the positional or/and orientational order is generally 
altered. In a simple picture, the full consistency of both potentials (if any) would be found near the bottom of the 
trap where the realistic potential becomes isotropic, i.e., for l/d^> 1. Since small A r favor grape configurations that 
violate both the positional and orientational orders of the parabolic case, the consistency should be rather searched in 
the opposite limit of strong A r ^> 1, where each trap should host exactly 3 colloids. In particular, the most promising 
phases are J- and A phases since they also obey the symmetry constraints of Vf eal . To avoid bending of extended 
trimers in these two phases, we also anticipate small screening, i.e., d « 1, otherwise a banana B shape might 
become more favorable. 
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FIG. 9: (Color online) Comparison of structures calculated within an isotropic parabolic potential (a) and realistic cosine 
potential Vf (b-d) with three neighboring magnitudes A r chosen to roughly correspond to the parabolic potential at the trap 
midpoint x/l — 0.5 (left panel). In (a) the parameter values are, from left to right: (l/d, kI, A, a) = (10,4,5,0.8) (J-* phase), 
(7, 2, 2, 0.8) (Q phase), (15, 4, 5, 1) (A phase), and (7, 4, 2, 1) (1Z phase). In (b-d) we use the same parameters as in (a) except 
Vf eal is used instead of parabolic trapping with A r — 20, 15, 10 for A = 5 and A r = 8, 6, 4 for A — 2. 




FIG. 10: (Color online) Single-trap calculations for parabolic (para) and realistic (real) trap potentials using A — A r = 50 
and a — 1 (left panel) and a = 0.8 (right panel). Shaded areas denote regions in parameter space where straight-trimer has 
lower energy then a banana B and grape Q clusters. The symbols are comparing full energy-minimization calculations for both 
trap potentials considering all neighboring traps: black disks denote exact consistency while the other symbols indicate that 
different phases are found with the 'realistic' and parabolic confinement potentials. 

Following these criteria and using a single-trap approximation (which is justified for small clusters), we identified 
regions in the (nlj/d) parameter space (see shaded areas in Fig. [TU| where one straight-trimer is energetically 
favorable against the banana B and grape Q structures. In Fig. [TU]we show such regions for realistic and parabolic 
potentials for a square and rectangular case (note that a has no effect in the parabolic approximation). The overlap 
between them gives an estimate where the consistency for J- (in a = 0.8) and A (in a = 1) should be searched within a 
proper calculation including all neighboring traps. We have checked this at few points of the phase diagram (symbols 
in Fig. [TU]) and indeed found consistency: denoted by black circles in Fig. [TO]are the results where the all-neighbour 



10 



calculation predicts the same phase (J 7 in Fig. [TOk and A in Fig. HOb ) for both potentials. These results agree with 
the single trap prediction. To the left from the overlapping region (at smaller l/d), the inclusion of the neighboring 
traps introduces bending of the trimers and thus prefers banana B phase over the J- or A (see green triangles in Fig. 
n~0|) . This effect is more visible for parabolic potential since it is stronger than realistic when A = A r . The above 
calculations were performed for a rather strong pinning strength (A = A r = 50); however, we have checked and found 
T and A consistency also for weaker trap potentials, with A = A r > 10. 

V. APPLYING AN EXTERNAL ELECTRIC FIELD 

Since most of the observed phases are generally close in energy, it is tempting to explore how an external perturbation 
could change the phase behaviour. We studied the response of the system to a static uniform electric field E with the 
in-plane direction. Our primary purpose was to numerically investigate the possibility of external pattern switching 
between the known phases, which might be appealing for experimental realizations. The ability to create and control 
colloidal crystal structures indeed has a wide range of applications in photonic and phononic materials, optical switches, 
photonic band-gap materials, and self-assembly of nanostructures 

In addition to either parabolic or cosine confinement, we use the following field potential at position (x,y): 

V E = ±Eq(xcosP + ysmfi)/l, (8) 

where Eq is a dimensionless field strength measuring the relative importance of external drag forces over Coulombic 
forces, and (3 is an angle between E and the longest lattice principal axis. The "±" sign selects between positively 
and negatively charged colloids, respectively. The amplitude Eq is considered small enough so that system acquires 
equilibrated ground state configuration. This requirement is easily fulfilled in a parabolic trap approximation where, 
in principle, we never reach stationary regimes with drifting colloids |2ll | even for largest values of Eq. However, the 
opposite may happen in a realistic trap confinement |2lJ where at critical Eq either cluster hopping or its breaking 
into smaller constituents takes place. We note that in the later case the critical field can be expressed analytically by 
E c ,2 = (1 + 2nd) exp(— 2nd) / (2d/l) 2 above which the two oppositely charged colloids (a dimer) unbind. The breaking 
of a trimer occurs at slightly lower fields that can be estimated numerically by E c ^ ~ 0.8-E Cj 2- 

We examined the influence of E on the selected phases shown in Sees. IIIII and IIVI using the EM and the gradient 
descent method. In Fig. [TT]we show the results (one-trap snapshots) for parabolic trapping. Apart from rather trivial 
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FIG. 11: Application of the external electric field E on colloidal clusters shapes in parabolic pinning regime for three field 
strengths Eo and directions depicted at the bottom of the figure. The other parameters are: (l/d, nl, A,a) — (10,4,5,0.8) in 
first row, (7, 4, 2, 1) in second row, and (7, 2, 2, 0.8) in third row. The asterisk symbol indicates that a cluster (when viewed 
periodically on a lattice) forms a percolated structure along the field direction. 

effects, such as bending (formation of the B shape out of the straight trimer, or a rocket 1Z deformation for Eq < 2), 
we identified also more interesting effects at strongest Eq = 5: a new shape evolved from 1Z as well as Q structure, 
and an elongation of the modified Q complex along the field direction lead to the formation of the percolated chain 
structure (denoted by asterisk symbol in Fig. Illl) . 

In a similar manner, we demonstrate in Fig. [T2]the influence of E on a particular A phase, again within a parabolic 
trap confinement; for Eq = 2 and an electric field angle (3 = or 7r/4, we identify a pattern switching to a (quasi) J 7 * 
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FIG. 12: Influence of an external electric field E on the A phase within the parabolic trap approximation for (l/d, kI, A, a) = 
(15, 4, 5, 1). The field directions and magnitudes are shown with arrows in the inset of each plot; arrow lengths correspond to 
magnitudes Eo — 0.5, 2, 5 of Eq. ([8]). For Eo — 2 and j3 = 0, 7r/4 (the middle column) we observe a pattern switching to the 
J-"* phase (actually, its close approximation, since the trimers here are slightly bent with ip\p ~ 0.9l7r). 



phase. Here, the transition A — > J-* appears to be on orientational basis, with more or less minor shape deformations. 
Since the trimers of field induced JF* phase are slightly shifted away from the trap center, the expected bending due 
to electric field is (almost) compensated by the restoring trap forces and the clusters appear almost straight. 

Wc have also explored the effect of the external field within the real trap confinement. As found by Reichhardt et al 
in Ref. [2l|, the colloids start to flow as soon as the electric field force becomes larger than the real trap confinement, 
but here we rather focused on static configurations at smaller field strengths. In fact, in a large enough system, any 
constant electric field would inevitably lead to ground states with uneven population of the traps, i.e. the trimers 
would prefer to move in the direction of the field to reduce the total energy of the system, leaving some of the traps at 
the opposite end empty. However, even though such states might have the lowest energy, they might not be accessible 
in the experiments and simulations. If, for instance, the system is first let to equilibrate at E = and then a small 
enough electric field is turned on, a local rather than global minimum is reached with colloids unable to surmount the 
large energy barriers on the way to the global minimum. This is the regime we are interested in here. 

To search for nearest local energy minimum with an increasing electric field, we introduced a gradient descent 
method (i.e., a method of steepest descent) on a 4-partite lattice (16 model parameters, see Fig. [T]) which brings a 
system continuously from an initial zero-field ground state to a finite-field metastable state using a finite step SEq 
(at fixed field direction j3). Depending on the amplitude of the increment SEq, the switching of the field may be 
interpreted as fast or slow, something that can also be tuned in an experimental setup and may affect the outcome. 
In our calculations we tried three step sizes, SE = 0.01,0.1,1, and found basically no differences in the resulted 
configurations. 

In Fig. [T3] we show the evolution of A and T phases with increasing E along two field directions using highest 
field resolution SE = 0.01. Focusing first on the A phase (upper row in Fig. Q2]), we identify at E ~ 3 (for both 
field directions) a rapid phase transition from a 2- to 1-partite phase, which can be also viewed as a pattern switching 
to an approximate J-* phase (with the internal angle being (p w 0.87T rather than 7r), similar to what we found in 
the parabolic confinement (middle column in Fig. I12[) . This tendency of reducing the "partiteness" of the phase by 
increasing the external uniform field seems natural. However, the opposite situation is observed in the evolution of 
the T phase, shown in Fig. ITBT c). where in a short interval, 3 < Eq < 4, the metastable phase adopts a higher, 
2-partite B structure. This rather surprising effect reflects the richness of the system studied, where the observed 
phases differ only slightly in energy and can therefore interchange easily under external perturbation. At larger field 
strengths again a uniform pattern takes place. 

We note that a similar field-evolution study in cosine trap confinement was made also for two representative Q and 
1Z phases with 6 colloids per trap, but we found no striking effects apart from trivial bending of the clusters, which 
was already observed in a parabolic confinement (see , e.g., Fig. 1111) . These results are included in a form of animated 
cartoons available as online supporting material [29(. 
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FIG. 13: (Color online) Local energy minima configurations represented by trimer angles &i (solid line) and ipi (dashed) as 
a function of external electric field Eo using the cosine confinement potential Eq. (0. The index i = 1...4 represents the 
four different values for a 4-partite lattice assumption. The upper row is calculated from initial (Eo = 0) A phase with 
(l/d, nl,A r /(2-K 2 ),a) = (15,4, 1.5, 1) for two field directions a) /3 = and b) /3 — 7r/4. In the lower row the initial (Eq — 0) 
state is the T phase with (l/d, kI, A r /(2iv 2 ), a) = (10,4,0.75,0.8) and field directions are again c) /3 — and d) f) — it/A. In 
all the plots we used field resolution SEo = 0.01. A few representative snapshots are added for illustration (see the supporting 
material [2S| for the animated cartoons). 



VI. CONCLUSION 



In summary, we have studied ordering of colloidal ionic trimers on square lattices. We observed a rich variety 
of crystal structures including crystals made of complex 6-colloid subunits. We have discussed their stability and 
the possibility for pattern switching by external control, which might be appealing for experimental realisations and 
applications. 

In case of the isotropic confinement potential, the symmetry is broken by the Yukawa interactions only. The 
orientational ordering is of 'intrinsic nature' and does not stem from an anisotropy present in the confinement potential 
itself. Three colloids can form a straight line or a bent banana shaped trimer. We have identified the parameter range 
where both are stable and the types of their orientational ordering. In the low n and small l/d part of the phase 
diagrams, where the inter-cluster interactions typically outweigh the confinement, the neighboring trimers cluster into 
larger complexes consisting of six colloids; their typical shapes being gra pe a nd rocket. At even smaller l/d and k, 
the colloids form percolated structures. Here, however, many-body effects [3fj] are expected to become significant and 
our assumption of the pairwise additivity of the colloid-colloid interaction becomes questionable [3ll . |32| ■ We have 
shown that most of the phases observed in the isotropic confinement persist also in the more realistic anisotropic 
external potential, provided that the confinement strength is large enough. However, their orientational ordering in 
the latter case is different due to the additional symmetry breaking by the laser field. The rich selection of the ground 
state structures and the relatively small energy differences between them enable efficient control over the structures 
by applying external fields. We have analysed the effect of a small external electric field and have shown that it can 
be used as a pattern switching tool. 

Not addressed here, but a subject for future work, is the thermal behaviour of the system. Another question lies 
is the effect of the interaction potential details on the stability of the observed structures, and in the many-body 
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interactions among colloids that arc expected to be important at high densities. Finally, in addition to systems with 
oppositely charged colloids, externally driven superparamagnetic colloids [33j and colloids interacting via Casimir 
interaction in critical binary mixtures [HI may offer instances where the present predictions -or adaptations thereof- 
could be realized. 

Acknowledgments: We acknowledge the support of the bilateral program Proteus of the Slovenian Research 
Agency and the French Ministere des Affaires Etrangeres et Europeennes, the support of the Slovenian research 
agency through the grant Pl-0055, and the support of the EU through the ARG (ERC-COLSTRUCTION 227758) 
and 7 th Framework Programme (ITN-COMPLOIDS 234810). 



[1] E.S. Hedges, Colloids, Hedges Press (2007). This monograph dates back to 1931 and offers an instructive overview of the 
state of the art at that time. 

[2] B. V. Derjaguin and L. Landau, Acta Physicochim. URSS 14 633 (1941); E. J. W. Verwey and J. T. G. Overbeek, Theory 

of the Stability of Lyophobic Colloids, Dover, Amsterdam (1997). 
[3] See e.g. D. Frenkel, Introduction to Colloidal Systems in Soft Condensed Matter Physics in Molecular and Cell Biology, 

edited by W.C.K. Poon and D. Andelman, Taylor and Francis (2006). 
[4] W.C.K. Poon, Science 304, 830 (2004). 

[5] M.E. Leunissen, C.G. Christova, A. P. Hynninen, CP. Royall, A.I. Campbell, A. Imhof, M. Dijkstra, R. van Roij and A. 

van Blaaderen, Nature 437, 235 (2005). 
[6] M. Brunner and C. Bechinger, Phys. Rev. Lett. 88, 248302 (2002). 

[7] C. Reichhardt and C.J. Olson, Phys. Rev. Lett. 88, 248301 (2002); C. J. Olson-Reichhardt and C. Reichhardt, J. Phys. 

A: Math. Gen. 36, 5841 (2003); C. Reichhardt and C. J. Olson-Reichhardt, Phys. Rev. E 71, 062403 (2005). 
[8] R. Agra, F. van Wijland and E. Trizac, Phys. Rev. Lett. 93, 018304 (2004). 

[9] A. Sarlah, T. Franosch and E. Frey, Phys. Rev. Lett. 95, 088302 (2005); A. Sarlah, E. Frey and T. Franosch, Phys. Rev. 
E 75 021402 (2007). 
[10] S. El Shawish, J. Dobnikar and E. Trizac, Soft Matter 4, 1491 (2008). 

[11] K. Mangold, P. Leiderer and C. Bechinger, Phys. Rev. Lett. 90 158302 (2003); J.-X. Lin, C. Reichhardt, Z. Nussinov, 
LP. Pryadko, and C.J.O. Reichhardt, Phys. Rev. E 73, 061401 (2006); J. Baumgartl, M. Zvyagolskaya, and C. Bechinger, 
Phys. Rev. Lett. 99, 205503 (2007); D. Deb and H.H. von Grunberg, J. Phys.: Condens. Matter 20, 245104 (2008); 
M. Schmiedeberg and H. Stark, Phys. Rev. Lett. 101, 218302 (2008). 

[12] E. Trizac, S. El Shawish and J. Dobnikar, Anais da Academia Brasileira de Ciencias 82, 87 (2010). 

[13] As will be discussed later, our model differs in some aspects from the one considered in Ref. [9J, where only discrete 
orientations were considered for the colloidal n-mers. Within such a framework, some approximations, backed up by Monte 
Carlo simulations, have been put forward to unravel the effect of temperature [9J. 

[14] M. Brunner and C. Bechinger, Progress in Colloid and Polymer Science 123, 156 (2004). 

[15] J.-P. Hansen and H. Lowen, Annu. Rev. Phys. Chem. 51, 209 (2000). 

[16] Y. Levin, Rep. Prog. Phys. 65, 1577 (2002). 

[17] E. Trizac, L. Bocquet and M. Aubouy, Phys. Rev. Lett. 89, 248301 (2002). 

[18] L. Belloni, Colloids and Surfaces A: Physicochem. Eng. Aspects 140, 227 (1998); J. Phys.: Condens. Matter 12, R549 
(2000). 

[19] It is worth stressing that the ground states we consider here are reached at constant k while the latter is a temperature 
dependent quantity. Such a point of view -consistent with the numerical approach of 0]- requires from an experimental 
perspective to consider large pinning amplitudes (strong laser fields) at constant temperature and constant electrolyte 
density. However, compared to the like-charge colloidal cases studied previously [6H11I]. there is no need here to tune 
correlatively the colloidal charge, since within a composite molecule in a given trap, the colloids are at contact: in the like- 
charge case, the extension of the molecules critically depends on the confinement amplitude, whereas here, this amplitude 
only affects the shape of the molecules, but not their extension. 

[20] J. Groenewold, W. Kegel, J. Phys. Chem. B. 105, 11702 (2001); A. Stradner, H. Sedgwick, F. Cardinaux, W.C. Poon, S.U. 
Egelhaaf, P. Schurtenberger, Nature 432, 492 (2004); A.I. Campbell, V.J. Anderson, J.S. van Duijneveldt, P. Bartlett, 
Phys. Rev. Lett., 94, 208301 (2005); S. Mossa, F. Sciortino, P. Tartaglia, E. Zaccarelli, Langmuir 20, 10756 (2004); S. 
Jiang, Q. Chen, M. Tripathy, E. Luitjen, K.S. Schweizer, S. Granick, Adv. Matt. 22, 1060 (2010). 

[21] C. Reichhardt and C.J. Olson Reichhardt. larXiv:0812.2447l 

[22] See e.g. Numerical Recipes, The Art of ScientiSc Computing, Cambridge University Press (1992). 

[23] We stress that our system exhibits at least two energy scales: the larger one is defined by the attractive Coloumb interaction 
between two touching colloids and is of the order l/2d, while the smaller one is the Yukawa interaction between the traps 
of the order exp(— nV)/l and is responsible for orientational order. The 4p-dimensional energy surface is therefore composed 
of deep minima each of which has again many shallow minima. 

[24] This mechanism somehow resembles the strain-induced colloidal clustering on periodic (triangular) substrates with more 
traps than colloids [35|]. Another somewhat related phenomenon is the appearance of the so-called cluster crystals in systems 
with soft repulsive potentials [36T439I] . although the present system is more complex and a delicate interplay between the 



14 



external confinement and the colloidal interactions governs its behaviour. 
[25] E. Trizac, L. Bocquet, R. Agra, J.-J. Weis and M. Aubouy, J. Phys.: Condens. Matt. 14, 9339 (2002). 
[26] R. Agra, E. Trizac and L. Bocquet, Eur. Phys. Journal E 15, 345 (2004). 

[27] At smaller l/d the checkerboard structure eventually becomes less favorable against a 2-partite structure: the corresponding 
energies, for instance, change from E 4p /K = -1.718, E 2p /K = -1.717 for {l/d, id) = (10,2) to E 4p /K = -0.9199, 
E 2p /K = -0.9203 for (l/d, kI) = (7,2). 

[28] A. Walther and A.H.E. Miiller, Soft Matter 4, 663 (2008). 

[29] See EPAPS Document No. []• 

[30] M. Brunner, J. Dobnikar, H. H. von Griinberg, and C. Bechinger, Phys. Rev. Lett. 92, 078301 (2004); 

J. Dobnikar, M. Brunner, H. H. von Griinberg, and C. Bechinger, Phys. Rev. E 69, 031402 (2004). 
[31] J. Dobnikar, Y. Chen, R. Rzehak, and H. H. von Griinberg, J. Chem. Phys. 119 (9) 4971 (2003). 

[32] D. Reinke, H. Stark, H. H. von Griinberg, A. B. Schofield, G. Maret, and U. Gasser, Phys. Rev. Lett. 98, 038301 (2007) 
[33] N. Osterman, I. Poberaj, J. Dobnikar, D. Frenkel, P. Ziherl, and D. Babic, Phys. Rev. Lett. 103, 228301 (2009). 
[34] C. Hertlein, L. Helden, A. Gambassi, S. Dietrich, and C. Bechinger, Nature 451, 172 (2008); 

F. Soyka, O. Zvyagolskaya, C. Hertlein, L. Helden, and C. Bechinger, Phys. Rev. Lett. 101, 208301 (2008). 
[35] S. Bleil, H. H. von Griinberg, J. Dobnikar, R. Castaneda-Priego and C. Bechinger, Europhys. Lett., 73 (3), 450 (2006). 
[36] B.M. Mladek, P. Charbonneau and D. Frenkel, Phys. Rev. Lett. 99 (23), 235702 (2007). 
[37] S. van Teeffelen, A. J. Moreno and C. N. Likos, Soft Matter 5, 1024 (2009). 

[38] N. Osterman, D. Babic, I. Poberaj, J. Dobnikar and P. Ziherl, Phys. Rev. Lett. 99, 248301 (2007). 
[39] J. Dobnikar, J. Fornleitner and G. Kahl, J. Phys.: Condens. Matter 20, 494220 (2008). 



